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^ Abstract 

We consider within QCD collinear factorization the process p + p — > jet + jet + X, 
where two forward high-pr jets are produced with a large separation in rapidity Ay 
(Mueller-Navelet jets). In this case the (calculable) hard part of the reaction receives 
CNj ■ large higher-order corrections ~ a™(Ay) n , which can be accounted for in the BFKL 

approach. In particular, we calculate in the next-to-leading order the impact factor 
(vertex) for the production of a forward high-p^ jet, in the approximation of small aper- 
ture of the jet cone in the pseudorapidity-azimuthal angle plane. The final expression 
for the vertex turns out to be simple and easy to implement in numerical calculations. 
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1 Introduction 



The production of two forward high-p^ jets in the fragmentation region of two colliding 
hadrons at high energies, the so called Mueller-Navelet jets [TJ, is considered an important 
process for the manifestation of the BFKL [2] dynamics at hadron colliders, such as Tevatron 
and LHC. 

The theoretical investigation of this process implies a combined use of collinear and BFKL 
factorization: the process is started by two hadrons each emitting one parton, according to 
its parton distribution function (PDF), which obeys the standard DGLAP evolution [3J. On 
the other side, at large squared center of mass energy i/i, i.e. when the rapidity gap between 
the two produced jets is large, the BFKL resummation comes into play, since large logarithms 
of the energy compensate the small QCD coupling and must be resummed to all orders of 
perturbation theory. 

The BFKL approach provides a general framework for this resummation in the leading 
logarithmic approximation (LLA), which means resummation of all terms (a s ln(s)) n , and 
in the next-to- leading logarithmic approximation (NLA), which means resummation of all 
terms a s (a s ln(s)) n . Such resummation is process-independent and is encoded in the Green's 
function for the interaction of two Reggeized gluons. The Green's function is determined 
through the BFKL equation, which is an iterative integral equation, whose kernel is known at 
the next-to-leading order (NLO) both for forward scattering (i.e. for t = and color singlet 
in the t-channel) [U E] and for any fixed (not growing with energy) momentum transfer t and 
any possible two-gluon color state in the t-channel [6]. 

The pro cess- dependent part of the information needed for constructing the cross section 
for the production of Mueller-Navelet jets is contained in the impact factors for the transition 
from the colliding parton to the forward jet (the so called "jet vertex"). 

Such impact factors were calculated with NLO accuracy in [7J, where a careful analysis 
was performed, based on the separation of the various rapidity regions and on the isolation 
of the collinear divergences to be adsorbed in the renormalization of the PDFs. The results 
of [7J were then used in [8] for a numerical estimation in the NLA of the cross section for 
Mueller-Navelet jets at LHC and for the analysis of the azimuthal correlation of the produced 
jets. This numerical analysis followed previous ones [9J [10] based on the inclusion of NLO 
effects only in the Green's functions. Recently we performed a new calculation [TT] of the jet 
impact factor, confirming the results of (TJ. 

In this paper we recalculate the NLO impact factor for the production of forward jets in 
the "small-cone" approximation (SCA) [T2jlT3]. i.e. for small jet cone aperture in the rapidity- 
azimuthal angle plane. Our starting point are the totally inclusive NLO parton impact factors 
calculated in [14] , according to the general definition in the BFKL approach given in Ref. |15j . 
The calculation is lengthy, but straightforward, since the standard BFKL definition of impact 
factor provides the route to be followed. The use of the SCA, moreover, allows to get a 
simple analytic result for the jet vertices, easily implementable in numerical calculations and 
therefore particularly suitable for a semi-analytical cross-check of the numerical approaches 
which treat the cone size exactly. 
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Figure 1: Diagrammatic representation of the forward parton impact factor. 



The paper is organized as follows. In the next Section we will present the factorization 
structure of the cross section, recall the definition of BFKL impact factor and discuss the 
treatment of the divergences arising in the calculation; in Section 3 we describe the procedure 
for the jet definition and the SCA; in Section 4 and 5 we present the details of the calculation 
at LO and NLO, respectively; in Section 6 we draw some conclusions. 



2 General framework 



We consider the process 

p(Pi) + P(P2) -> jet(fci) + jet(fc 2 ) + X (1) 

in the kinematical region where the jets have large transverse moment 80, k 2 ~ k 2 ^> Aq CD . 

This provides the hard scale, Q 2 ~ k 2 2 , which makes perturbative QCD methods applicable. 
Moreover, the energy of the proton collision is assumed to be much bigger than the hard scale, 
s = 2pi ■ p 2 > k 2 2 . 

We consider the leading behavior in the 1/Q-expansion (leading twist approximation). 
With this accuracy one can neglect the masses of initial protons. The state of the jets can be 
described completely by their (pseudo)rapiditie9^]?/i i 2 and transverse momenta k\ j2 . Moreover, 
we denote the azimuthal angles of the jets as <pi j2 . 

In QCD collinear factorization the cross section of the process reads 

i i 

— — = f I dxidx2fi(xi, /ip)fj(x 2 , LXp) ^( XlX2S J t-E- ; (2) 

dy 1 dy 2 d 2 k 1 d 2 k 2 ij=g J J ' dy 1 dy 2 d 2 k 1 d 2 k 2 

where the i,j indices specify parton types, i,j = q,q,g, fi(x,jj,p) are the proton PDFs, 
the longitudinal fractions of the partons involved in the hard subprocess are x\ t2 , fxp is the 
factorization scale and da(xix 2 s, /xf) is the partonic cross section for the jet production. 



1 See Eq. ((3]) below for the definition of the transverse part of a 4- vector. 

2 For massless particle the rapidity coincides with pseudorapidity, y — rj, the latter being related to the 



particle polar scattering angle by r\ — — In tan 



2 ■ 
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Figure 2: Parton-Reggeon collision, the jet is formed by a single parton. 



It is convenient to define the Sudakov decomposition for the jet momenta, 



k 2 

h = aipi + — p 2 + h,± , k 2 ± = -k 2 , (3) 



k 2 



-Pi + OL2V2 + k 2 ,± , k 2 ± = —k 



_. _ 2 > 

where the jet longitudinal fractions «i 5 2 are related to the jet rapidities by 

1 a\s 1, a^s . 

Sl = y 2 = - 5 .n ¥ , (4) 

and = d?/2 = — ^ i n the center of mass system. 

We consider the kinematics when the interval of rapidity between the two jets, 

Ay = yi - y 2 = In , (5) 

is large. Since the jet longitudinal fractions are equal or smaller (in the case of additional 
QCD radiation) than the ones of the participating partons, ai < x\, a 2 < x 2 , we are in 
a situation where the energy of the partonic subprocess is much larger than jet transverse 
momenta, Xix 2 s ^> k 2 2 (k 2 and k 2 2 are considered to be of similar order ~ k 2 ). In this 
region the perturbative partonic cross section receives at higher orders large contributions 
~ a™ ln n J^, related with large energy logarithms. It is the aim of this paper to elaborate the 
resummation of such enhanced contributions with NLA accuracy using the BFKL approach. 

Let us remind some generalities of the BFKL method. Due to the optical theorem, the 
cross section is related to the imaginary part of the forward proton-proton scattering ampli- 
tude, 

1m s A 

a = —— . (6) 
In the BFKL approach the kinematic limit s 3> k 2 of the forward amplitude may be presented 



3 



in D dimensions as follows: 



8+ioo 

s f dP~ 2 q f dP~ 2 q f duo ( s 

5— too 

where the Green's function obeys the BFKL equation 

w g 2 ) = 5 D " 2 (gi - q 2 ) + / d D ~ 2 qK(qi, q) G u (q, <fi) . 



(7) 



What remains to be calculated are the NLO impact factors $1 and $2 which describe the 
inclusive production of the two jets, with fixed transverse momenta ki, k 2 and rapidities y\, y 2 , 
in the fragmentation regions of the colliding protons with momenta p\ and p 2 , respectively. 
The energy scale parameter s is arbitrary, the amplitude, indeed, does not depend on its 
choice within NLA accuracy due to the properties of NLO impact factors to be discussed 
below. 

For definiteness, we will consider the case when the jet belongs to the fragmentation region 
of the proton with momentum pi, i.e. the jet is produced in the collision of the proton with 
momentum p\ off a Reggeon with incoming (transverse) momentum q and denote for shortness 
in what follows its transverse momentum and longitudinal fraction by k and a, respectively. 

Technically, this is done using as starting point the definition of inclusive parton impact 
factor, given in Ref. [H], for the cases of incoming quark(antiquark) and gluon, respectively 
(see Fig. 1). Here we review the important steps and give the formulae for the LO parton 
impact factors. 

Note that both the kernel of the equation for the BFKL Green's function and the parton 
impact factors can be expressed in terms of the gluon Regge trajectory, 

j(t) = 1 + u(t) , (9) 

and the effective vertices for the Reggeon-parton interaction. 

To be more specific, we will give below the formulae for the case of forward quark impact 
factor considered in D = 4 + 2e dimensions of dimensional regularization. We start with the 
LO, where the quark impact factors are given by 

^ 0) (?) = E / r 8?tf) [r £ 3 ^* dPa ' (10) 

{a} 

where q is the Reggeon transverse momentum, and T^j denotes the Reggeon-quark vertices in 
the LO or Born approximation. The sum {a} is over all intermediate states a which contribute 
to the q — >■ q transition. The phase space element dp a of a state a, consisting of particles with 
momenta £ n , is (p q is initial quark momentum) 

d Pa = (2,f ^ L + q-Y, In) If So\ • (U) 

\ n£a / nGa ^ ' n 
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while the remaining integration in (ITUj) is over the squared invariant mass of the state a, 

M 2 a = ( Pq + qf . 



In the LO the only intermediate state which contributes is a one-quark state, {a} = q. 
The integration in Eq. ( flOl) with the known Reggeon-quark vertices T q °q is trivial and the 
quark impact factor reads 

*Tti)=9 2 ^ 1 , (12) 
where g is QCD coupling, a s = g 2 /(4n), N = 3 is the number of QCD colors. 

In the NLO the expression (flOj) for the quark impact factor has to be changed in two ways. 
First one has to take into account the radiative corrections to the vertices, 

r(o) v r _ r(o) , r (i) 

1 qq ~ L qq L qq ~ L qq • 

Secondly, in the sum over {a} in (Till]) , we have to include more complicated states which 
appear in the next order of perturbative theory. For the quark impact factor this is a state 
with an additional gluon, a = qg. However, the integral over M 2 becomes divergent when 
an extra gluon appears in the final state. The divergence arises because the gluon may be 
emitted not only in the fragmentation region of initial quark, but also in the central rapidity 
region. The contribution of the central region must be subtracted from the impact factor, 
since it is to be assigned in the BFKL approach to the Green's function. Therefore the result 
for the forward quark impact factor reads 



4,(9 , s »> = (|IJ E J "af r ««(9 > iW? )1* - M ») 

4/^H^ (W "*K<^) ■ (13) 

The second term in the r.h.s. of Eq. ffT31) is the subtraction of the gluon emission in the 
central rapidity region. Note that, after this subtraction, the intermediate parameter s\ in 
the r.h.s. of Eq. ( fl3|) should be sent to infinity. The dependence on s\ vanishes because of the 
cancellation between the first and second terms. is the part of LO BFKL kernel related 
to real gluon production, 

The factor in Eq. ( !T3|) which involves the Regge trajectory arises from the change of energy 
scale (q 2 — > Sq) in the vertices T. The trajectory function oj(t) can be taken here in the 
one-loop approximation (t = —q 2 ), 

9 2 t N f d D ~ 2 k 2 , r r(i - £ ) r 2 (eK 2NF 
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In the Eqs. (jTOl) and (TT3]) we suppress for shortness the color indices (for the explicit form 
of the vertices see [H]). The gluon impact factor <& g (q) is defined similarly. In the gluon case 
only the single-gluon intermediate state contributes in the LO, a = g, which results in 

*?\?) =^ ] (th (16) 

here Ca = N and Cp = {N 2 — 1)/(2N). Whereas in NLO additional two-gluon, a = gg, and 
quark- ant iquark, a = qq, intermediate states have to be taken into account in the calculation 
of the gluon impact factor. 

The definition of inclusive parton impact factors involves the integration over all possible 
intermediate states appearing in the parton-Reggeon collision. Up to the next-to-leading 
order, this means that we can have one or two partons in the intermediate state. Then, 
in order to allow for the inclusive production of a jet, these integrations must be suitably 
constrained to take into account that the kinematics of the parton or the pair of partons 
which generate the jet is fixed by the jet kinematics. 



3 Jet definition and small-cone approximation 

At LO the (totally inclusive) parton impact factor takes contribution from a one-particle 
intermediate state; equivalently, only one parton is produced in the collision between the 
incoming parton and the Reggeon, as shown in Fig. [2j Therefore, the kinematics of the 
produced parton is totally fixed by the jet kinematics. At NLO we have both the virtual 
corrections (which have the kinematical structure shown in Fig. EJ) and also two-particle 
production in the parton-Reggeon collision. The jet in the latter case can be either produced 
by one of the two partons or by both together. If we call the produced partons a and b, we 
have the following contributions, as shown in Fig. [3] (see, for instance, Ref. [16]): 

1. the parton a generates the jet, while the parton b can have arbitrary kinematics, provided 
that it lies outside the jet cone; 

2. similarly with a -H- b; 

3. the two partons a and b both generate the jet. 

The cases 1. and 2. are replaced in the actual calculation by the following two (as 
illustrated in Fig. 0J: 

1. the parton a generates the jet, while the parton b can have arbitrary kinematics ("in- 
clusive" jet production by the parton a); then, the case when the parton b lies inside 
the jet cone is subtracted; 

2. similarly with a -H- b. 
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Figure 3: Parton-Reggeon collision, two partons are produced and the jet is formed either by 
one of the partons or by both partons. 



Let us introduce now the "small-cone" approximation (SCA). In view of the discussion 
above, we should define it in the two cases of jet generated by one parton or by two partons. 

The relative rapidity and azimuthal angle between the two partons are 

q-k-k 2 



2 (2 k 2 



A0 



arccos • 



\k\ \q — k\ 



C = i - C - 



Let the parton with momentum k and longitudinal fraction ( generate the jet, whereas the 
other parton (with momentum q — k and longitudinal fraction () is a spectator. We introduce 
the vector A such that 

Then, for A — > we have 

^ = s l^-^l), Ay «*"■*> 




fc 4 J C k 2 

thus the condition of cone with aperture smaller than R in the rapidity-azimuthal angle plane 
becomes 

C 2 A 2 

A0 2 + A y 2 = L^<^ 
C k 2 

and therefore 

|A| < ^\k\R. 

The situation is different when both partons form a jet. In this case the jet momentum 
is k = k\ + k 2 and the jet fraction is 1 = ( + (. The relative rapidity and azimuthal angle 
between the jet and the first (second) parton are 

1 , k 2 Ai k-h 

Awi = — m — — , A0i = arccos — — — , 

2 (2f,2 Vl \ kl \\ k \ 
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a=jet 



jet 





Figure 4: The production of the jet by one parton when the second one is outside the cone 
can be seen as the "inclusive" production minus the contribution when the second parton is 
inside the cone. 



Ay 2 = - In v — ' 

2 -(2 k 2 



A0 2 



arccos • 



k- (k-ki) 
\k\\k — ki\ 



Introducing now the vector A as 
we find 

Ay 2 + A0 2 



fci = (k + A 
A 2 



Ay\ + A<& 



( 2 k 2 C 2 k 2 
so that the requirement that both partons are inside the cone is now 



|A| < R\k\ min(CC) 



4 Impact Factor in the LO 

The inclusive LO impact factor of proton may be thought of as the convolution of quark and 
gluon impact factors, given in Eqs. ( I12|[T()|) . with the corresponding proton PDFs, 

= C dx (^f g (x) + fa(x)) , C = 9 2 ^f^ = 2t«.7^ ■ (17) 

\ a=g,g / 

In order to establish the proper normalization for the jet impact factor, we insert into the 
inclusive impact factor ( fl~T|) the delta functions which depend on the jet variables, transverse 
momentum k and longitudinal fraction a: 




S 



In what follows we will calculate the projection of the impact factor on the eigenfunctions 
of LO BFKL kernel, i.e. the impact factor in the so called (is, n)-representation, 

^n) = / d 2 q^^= (qT~" ■ (19) 
J q 7TV 2 

Here is the azimuthal angle of the vector q counted from some fixed direction in the trans- 
verse space. 



5 NLO calculation 



We will work in D = 4 + 2e dimensions and calculate the NLO impact factor directly in the 
(is, n)-representation f|T9|) . working out separately virtual corrections and real emissions. To 
this purpose we introduce the "continuation" of the LO BFKL eigenfunctions to non-integer 
dimensions, 

{qye^^{qy^(q.l) n , (20) 

where 7 = iis — ■= and 1 2 = 0. It is assumed that the vector I lies only in the first two of 
the 2 + 2e transverse space dimensions, i.e. / = e\ + ie 2 , with e 2 2 = 1, e*i • e 2 = 0. In the 
limit e — > the r.h.s. of Eq. (|20|) reduces to the LO BFKL eigenfunction. This technique was 
used recently in Ref. [T7]. An even more general method, based on an expansion in traceless 
products, was uses earlier in Ref. [T5] for the calculation of NLO BFKL kernel eigenvalues. 
In the case of interest, I 2 = 0, these two approaches lead, actually, to similar formulas. 

Thus, for the case of non-integer dimension the LO result for the impact factor reads 

m j c p+2e k / c 



q 2 k 



C da 



2 



5^(k-q)l^f g (a)+J2fa(oc)) , (21) 



a=q,q 



which in the (v, n)-representation gives the result 
ttV2 k 2 d& (is, n) (C A 




/*(")+£/■(«) (* 2 ) 2 {k-l) ■ (22) 



C dad 2 ^k \C F - 

Collinear singularities which appear in the NLO calculation are removed by the renormal- 
ization of PDFs. The relations between the bare and renormalized quantities are 



f q (x) = f q (x, fi F ) - a (| + In J f [P qq (z)f q (% fi F ) + P qg (z)f g (%n F ) 

X 

f g (x) = f g (x,n F ) - f (| + lnf ) f± [P„(z)f 9 (*,H F ) +P 99 (z)f 9 (%»F)} , (23) 
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where 



| = ~ + 7£ — ln(47r) f=a ^fe^ , and the DGLAP kernels are given by 



T R [z 2 +\l - z) 2 ] 



1C A 



1 - z 
1 



1 + z 2 



-6(1 - z) 



1 



+ 



6(1 -z) 



(24) 
(25) 

(26) 
(27) 



with Tr = 1/2. Here and below we always adopt the MS scheme. 



Now we can calculate the collinear counterterms which appear due to the renormalization 
of the bare PDFs. Inserting the expressions given in Eqs. (123]) into the LO impact factor (1221) . 
we obtain 



Try/2 k 2 d® J (v, n) \ co \\ lne&v c , t , _ _« ! / 1 
C dad 2+2 4 ~ 2tt \e 



+ In 



/4 



k-l 



dz 

z 



(28 



E (p-w/. (f ) + PMf. (§)) + § (^ma (f ) + E /• (f )) 

a=g,<j \ a=q,q / 



The other counterterm is related with the QCD charge renormalization, 

11C A 2n f 



a s = a s (/jL R 
and is given by 



a a (p R ) 



^ + ln4 



Bo 



(29) 



7T 



y/2Pd$ J (v 



> ^7 (charge c.t. *-^s 



dad 2+2 'k 



2tt V e 



^ ~ + ln^ fc 



/4 



x 



\a=q,q 



Up 



fc-Z 



n F d? 

' —6(1 -z) 



11C A _ n f 

~~6 y 



(30) 



To simplify formulae, from now on we put the arbitrary scale of dimensional regularization 
equal to the unity, fi = 1. 

In what follows we will present intermediate results always for 7rv ^^ 2 d ^ .} U '™1 , which we 
denote for shortness as 

nV2k 2 d<S> J (v,n) = 
C dad 2+2 4 ~ 

Moreover, a s with no argument can always be understood as a s (jj>r)- 



(31) 



We will consider separately the subprocesses initiated by a quark and a gluon PDF, and 
denote 

I = I q + I 9 . (32) 
We start with the case of incoming quark. 
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5.1 Incoming quark 



We distinguish virtual corrections and real emission contributions, 

I q = I v q + J* . (33) 

Virtual corrections are the same as in the case of the inclusive quark impact factor, there- 
fore we have 

q 2tt (47r) e eT(l + 2e)l / V J J ( v s; Ja \( J 



a=q,q 



X { CF {l-TTTe +1 )- "' (1 + ^3 + 20 + C * P + - f) - 2 *' M + 

+ 4(l + 2e)(3 + 2e) ~~ e(l + 2e) ~ 4(1 + 2e) ~ 2) } ' ^ 

/ W-<? 2 ) 

Note that the contribution ~ In || in Eq. ( 134|) originates from the factor ( |l ) in the 

definition of the NLO impact factor, see Eq. (Tl3"j) . which is accounted for virtual corrections 
in the BFKL approach. 

We expand ( 1341 in e and present the result as a sum of the singular and the finite parts. 
The singular contribution reads 



a 



a=q,q 



2 o \ n f , ^< s o , 11 



whereas for the regular part we obtain 



1 



5n/ „ /85 vr 1 



x ^ 8Cf + T" Ci U + T U • (36) 



Note that + differs from 1^ by terms which are 0(e). 



5.1.1 Quark-gluon intermediate state 

The starting point here is the quark-gluon intermediate state contribution to the inclusive 
quark impact factor, 

*™ - *sr 0% M> -So- + CA P-^® ' (37) 
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where Pi and 02 are the relative longitudinal momenta (fl\ + 02 — 1) an d k\ and &2 are the 
transverse momenta (ki + fc 2 = <z) of the produced gluon and quark, respectively. 

We need to consider separately the "inclusive" situations when either the quark or the 
gluon generate the jet, with the kinematics of the other parton taken arbitrary. We denote 
the corresponding contributions as I^ q and I^ g , 



rR 



jR + p 



rR 



We start with the case of inclusive jet generation by the gluon, 1^ . 
a) gluon "inclusive" jet generation 

The jet variables are k — k\, ( = 0\ (02 = ( = 1 — (, k 2 = q — k), therefore we have 



q;g 



27r(4vr) 



d 2+ 



2c 



71 



1+e 



q-l 



^ V f ( - 

a=q,q ^ 



l + c 2 + < 2 



Q- 



+ C A 



^r-k-q 



(q- k) 2 (q- 



C 



(3* 



It is worth stressing the difference between the previous calculations of NLO inclusive 
parton impact factors and the present case of production of a jet with fixed momentum. 
In the parton impact factor case, one keeps fixed the Reggeon transverse momentum q and 
integrates over the allowed phase space of the produced partons, i.e. the integration is of the 
form J 2((i-() d 2+2e k . . . In the jet production case, instead, we keep fixed the momentum of the 

parent parton (, k, and allow the Reggeon momentum q to vary. Indeed, the expression (|38p 
contains the explicit integration over the momentum q with the LO BFKL eigenfunctions, 
which is needed in order to obtain the impact factor in the (is, n)-representation. 

The q- integration in fl38|) generates 1/e poles due to the integrand singularities at q — > 
k/( for the contribution proportional to Cp and at q — > k for the one proportional to Ca- 
Accordingly we split the result of the ^-integration into the sum of two terms: "singular" and 
"non- singular" parts. The non-singular part is defined as 



a s f <K_ 
2rt(ArtyJ C 



a=q,q 



C A 



c v c 



d 2+2 *q tf-k-q 



7T 



l+e 



k) 2 [q-l 



■2-(l7TY J <, ^J a {c)° A C 



a=q,q 



a 



c fi + e + ee 



c 



k-l 



k 2 ) 7+e ~* Ik- 1 
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where n is a unit vector, n 2 = 1. Taking this expression for e = we have 



(I R ) 



a. 



where we define the function 



k ■ I 




a=q,q 



c; cv c 



a, 



h = Ji(n,7,C) 



71 (a-n) 2 (a-f)' 

with the azimuth of the vector a counted from the direction of the unit vector n. 
For the singular contribution we obtain 



a j[l- e ] i r 2 (l + e) ft.2V +e 
2tt (47r) e e r(l + 2e) 



k 



k-l 




7 



a=q,q 



l + c 2 + < 2 



c 



r(l + 2e)r(f- 7 -e)r(f + l + 7 + e) 



cv „: v ~ ;:;r v2 ,„l v 1 ' ^, c 2£ ~ 27 + c A 



r(l + e)r(l-e)r(f- 7 )r(f + l + 7 + 2e) 



Expanding it in e we get 



(I R ) 



a,r[l-e]r 2 (l + e) ^ 2 ^T+ e -§ 



2tt e(47r) e T(l + 2e) 

x (^(0 



k-l 




dC 

c 



+e 



i + C 2 



^ + c 27 



^r 27 (x(n,7)-21nC) + 2C A ln^ 



(39) 



+ c(cvr 27 + c A ) 



where 
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(40) 



/ 71 \ { 71 

X(n, 7 ) = 2^(1)-^- -7j -^(2 +1 + 
is the eigenvalue of the LO BFKL kernel, up to the factor Na s /ir. 
b) quark "inclusive" jet generation 

Now the jet variables are k = k 2 , ( = P>2 (Pi = C> ^1 = Q~k). The corresponding 
contribution reads 



I?. 



R 



m 27r(47r) e J 7r 1+e 



d 2+2e <f _ 



(9 2 ) 1_t (?• f )" / f £ /. 
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X 



i + C 2 + < s 
(i-C) 



+ c A 



(41) 



We will consider separately the contributions proportional to Cp and Ca- 
bi) quark "inclusive" jet generation: CV-term 

Note that the integrand of the CV-term is not singular at £ — > 1. We use the decomposition 

/ 



k' 



k- 



(q-k) 2 (q-l) (q-k) 2 +(q-l) 



1 1 

+ 



{q-kf 



in order to separate the regular and singular contributions. The regular part is given by 

i 



^ ( P r ^-, T /f E/ 

a 0=5,5 



«\ C A + C J + eC 



c 



xC, 



d 2+2 'd 



l+e / ,\2 

(a-n) 2 + (a-f 



♦2\7-q / a- J 



a — 71 



- 1 (a 
+ — 



2^7-f ( 



-27 



a- 



n\2 



Cf ®s ( ^2 



Therefore for e = we have 

W)r 2 7T 

where we define the function 

d 2 a 



k-l 




a=q,q 



«\ ca + c 2 



c; c 2 



h = I 2 (n,7,C) = 



7T 



[a — n 



) 2 + (*-?)' 



(a 2 )V n <^-l , (a 2 )V n *-C~ 
(a — n) ; 



27 



+ 



(«-f) 2 



(42) 



(43) 



(44) 



The singular part is proportional to the integral 



/ 



d 2+2 *q 



k' 



-l+t -, / -\ 2 

(? - A;) 2 + (q 



C 



(q- kf 



(q- W 



r(i- e )r 2 (i + e ) fr2 \T+*-Sf r _ r WC 2 



e r[i + 2e] v ; v ■ ) VC 

therefore, for the singular part of the CV-term we have 



a s r(l-e)r 2 (l + e) f d( , p 



2tt e(4vr) e T(l + 2e) J ( 



7+e-i 



a=q,q 



a 
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xCf 



;i+c 2 +< 2 ) fc 



(0 26 ( 1 + ^ 27 ) 



(i-C) 

The next step is to introduce the plus-prescription, which is defined as 

11 a 



no 

(i-0- 



(i-0 



(i-0 



(45) 



for any function F((), regular at ( = 1. Note that 



(i-0 



2e-l 



(i-or 1 + ^(i-o = 7r^(i-o + 



2e 



2e 



(i-0- 



+ 2e 



Mi - 
i-C 



+ 0(t 



Using this result, one can write 

(i + C 2 + < 2 ) fC 



c f 



(i-0 



(1 + r 27 ) = c F 



^d-o + ^(i + c-) 



n;i r 27 )(c+2(i + c 2 )f i 4 L ^l -2(i + c 2 ^ lnC 1 -" 2i 



v i-c 



(i-0 



+ 0(e 2 



\ - 3) 5(1 - C) + + ^(1 - 0) (1 + r 27 ) + 0(e) 

= C F (- e - 3) 5(1 - C) + P„(C) (1 + C 27 ) + 0(e) . 

Taking this into account and expanding in e the singular part of the (Tp-term, one gets the 
following result for the divergent contribution: 



a=q,q 



(t r )C F = q s r[i-6] r 2 (i + e) ^ 

\ m) s 2vr e(4irY T(l + 2e) V 

x \c F - 3) 5(1 - c) + P„(C) (1 + C 27 ) 

+ e ^(i + r 27 ) fc + 2(i + c 2 )f^ i ^'l -2(i + c 2 ^ lnC 



i-C 



(i-O 



(46) 



b 2 ) quark "inclusive" jet generation: CVterm 

The C^-contribution needs a special treatment due to the behavior of (l4Tj) in the region 
C — > 1. We use the following decomposition: 



c A 




_ C A 2 1 


(i-O 
h Ca \ 


\(q-ky(q-iy J 
1 1 


2 (g - (1 - C) 


2(1-0 


(9-f) 2 (q-ky \C] ^-%y(q-\f 
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The second term in the r.h.s. is regular for ( — > 1 and can be treated similarly to what we 
did above in the case of the CV-contribution. The first term is singular and the integration 
over ( has to be restricted, according to definition of NLO impact factor, see Eq. (|T3l) . by the 
requirement 



M 2 QG < s A 



M 2 QG 



k 2 k 2 

_L + _! _ ff 2 



(g - ky 

i-C 



and assuming the sa parameter to be much larger than any scale involved, sa ^> q 2 ,ki 2 - 
Therefore the C integral has the form 

i-Co 



no 

i-C 



for Co 



(g - 

SA 



-»■ 



Using the plus-prescription ( 1451) one can write 

i-Co 



no 

i-C 



i-0- 



F(l)ln- 

{,0 



for Co -> 



(47) 



for any function -F(C) riot singular in the limit £—>■!, and 







\ 



(q-k] 



ln 



A;) 2 (q-ky 



c A 



2 (g*_ife)2(l-0 + 2(1-C) 



{q-k) 2 {q-±) 2 



We remind that the definition of NLO impact factor requires the subtraction of the contribu- 
tion coming from the gluon emission in the central rapidity region, given by the last term in 
Eq. (1131) . which we call below "BFKL subtraction term" . After this subtraction the parameter 
sa should be sent to infinity, sa — > oo. Our simple treatment of the invariant mass constraint, 
Mq G < sa, anticipates this limit sa — > oo, therefore we neglect all contributions which are 
suppressed by powers of 1/sa- Moreover, the first term in the r.h.s. of the above equation 
should be naturally combined with the BFKL subtraction term, giving finally 







q 2 -k-q^ + k 



\ (q-k) 2 [q 



-3*(i-C)-^- 

2 V ^(q-k) 2 



In 



s 



(g- k)' 



c 



2 tf- k ) 2 (1 - C)+ 2(1-C) 



(g- 



k\2 



(g - ky 



C 2 

c 



(q-k) 2 (q-^) 2 



a result where the artificial parameter sa cancels out, as expected. 
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After that, we are ready to perform the q- integration, which naturally introduces the 
separation into singular and non-singular contributions. The singular contribution reads 



o.r ( i- e )r>(i +e ) rpr~*(i.TY ff^f-Cz 

a=q,q ^ 



2tt e(4vr) e T(l + 2e) V / 



X 



x CU h , ,2 , - 2 , f r(l + 2e)r(f-7-e)r(f + l+ 7 + 6) 
2 I +^ +e ^\r(l + e)r(l-e)r(f-7)r(| + l + 7 + 26) 

5(1-0 (lng + ^(|-7-e)+^(l + 7 + | + 2e)-^( e )-^(l)) - 
(c -2e- 27 _ 1) - 



(i-0- 



+ 



i-C 



Expanding this expression in e and using that 

(C - 2£ -2 7 _ 1} (C -2e-2 7 _ 1} 



i-C 



and 



^2e-l ^-2e + £-2 7 -2ej = ^-2e + £-2 7 -2e^ ^ 



(i-C) + 
5(i-C) 



2e 



(i-0- 



0(e] 



we get the divergent term 



(i r ^ _ r(i - e) r 2 (i + e) ( p y«-$ (? ry [ <K r, f« 

l J W s - ^ e(4 ^ r(i + 2e) r J r y JT^ u \c 



x^C^(l-C)lng 



5(1 - ( X(n, 7) In g + - (V (l + 7 + |) - ^ " t) ~ xV, 7 ) 

ln( 



+ (i + c 2 )((i + r 2 ')(^| 7 -(i^ 



+ 



(i-0 



(4* 



The regular contribution differs from (142 p only by one factor and reads 



V 2 7T V 



fc-z 



a=q,q 



a\ C fl + C 



(J C V c 



c) both quark and gluon generate the jet 

In this case the jet momentum is k = h\ + k 2 and the jet fraction is 1 = C + C- Introducing 
the vector A as 

k l = (k + A , 
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the contribution reads 



a 



I R = 
m+g 2 7t(4 



-(,-r(-TE/.( Q )/^/ 

a=q,q ' n 



X 



i + C 2 + < 2 



(k 2 ((k-A + A 2 ) 



A 2 (CA; + A) 2 A 2 ((k + A) 2 ((k - 



In the small-cone approximation (SCA) we need to consider only 



A 2 . 



a=q,q n 



7T 1+e A2 



where |A max | = \k\Rmin((, £). Using that 



A 2 



d 2+2e A 1 



1+e A 2 er(l + e) 



we get 



I R = 



a.r(l-e) T 2 (l + e) /~.,\'+<- 



q > q+9 2tt e(4vr) e T(l + 2e) 



er(l + 2e) 

-,\n — 
k-l) J2fa^)C F R 



2c 



a=q,q 



x / dC(min(C,C)) 



2e i + C 2 + < 2 



x 



c 



„ a a r(i-e) r 2 (i + e ) / p y+-f /r r - 

~ 2tt e(47r)^ r(l + 2e) ' ' 
Y,fa{a)C F R 2 ^ 



a=q,q 



7-2 + f (2-T + 3hl2 J 



d) gluon "inclusive" jet generation with the quark in the jet cone 

— * 

We introduce the vector A such that 

— * 

— * — * 

where k coincides with ki, the transverse momentum of the gluon generating the jet. 
contribution reads 



a. 



27r(47r)' 
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A 2 

x i + e + ,e ft^JL, (52) 

C J TT 1+t A 2 



where now |A max | = \k\R$. We get 



2< 



X /fc-E.(f)(f)^ 



1 



a s T(l-e) T 2 (l + e) 
'2tt e(47r) e T(l + 2e) 

x(p ot (C) l + 26ln^ 



„ ^ «=9:9 ^ 



Note the overall minus sign, which means that this contribution is a subtractive term to the 
gluon "inclusive" jet generation. 

e) quark "inclusive" jet generation with the gluon in the jet cone 

In this case we have 

with k identified with hi, the transverse momentum of the quark generating the jet. The 
contribution reads 



q;q,-g 



2tt(47 



a=q,q 



(i-C) 

where |A max | = We get 



A 2 

„ — „ max — * 

l + C 2 + < 2 /• ^ 2+2e A 1 

x 



f d A 1 

y — ^ ■ (54) 



fl m„ Til O l'-(l ■ Q ( , M/ rvv , ry2i 

q '' q - g 2tt e(47r) e T(l + 2e) 



i 

x 



/f-E/.(!)(9 W 

a. r(l - e) + e) ( - 2 pi (£ . f) » fl2< | ^ g 



2tt e(4vr) e IY1 + 2e) 
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5.1.2 Final result for the case of incoming quark 



Collecting all the contributions calculated in this Section and taking into account the PDFs' 
renormalization counterterm fl28|) and the charge counterterm fl30|) . we find that all singular 
contributions cancel and the result is 



Iq ~27T { ' 



k-l 



a=q,q 



(56) 



{P m (0 + ^PUO ) ^ K ~ 2C- 27 In R {P m (0 + P gq (0} - T ln ^ ~ 
+CaS(1 - C) [x(n, 7) In g + ^ + y + \ (V (l + 7 + |) - V Q - t) " xV 7)) 



+2 



(i + r 27 )x(n, 7 ) 

2(1-0+ 
ln(l - C) 



ln(l - C) 

i-C 



+C A ( + C F ( + 



i-C 

1 + C 2 f . 



+ <J(1-C) \ C F 31n2- 



7r 2 9\ 5n/ 



21n C 

.c 2 2 i-c 



3 2 



9 



c 



C A ±h + 2C A In ^ + C F C 2l (x(n, 7) - 2 In C) 



5.2 Incoming gluon 



We distinguish virtual corrections and real emission contributions, 

/ =I V + I R 



Virtual corrections are the same as in the case of inclusive gluon impact factor, 



jV = a s r[l-e] i r 2 (l + e) /pV 



2tt (Any eV(l + 2e)V" J \ k ' 1 ) 



C f 



11 + 9e 



£2 e 2(l + 2e)(3 + 2e) 



+ V»(l-e) -2^(l + e)+V(l) 



+ ;i + e)(l + 2e)(3 + 2 e )J + n/ ^ (1 + e)(l + 2e)(3 + 2e) J j ' 
Expanding it in e we obtain the following results for the singular, 

2vr (4vr) e e r(l + 2e)V / V / JbK ' C F 



(57) 
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k 2 e 6/ 3 



and the finite parts, 



(a-£(pr(j-o"«-)§h@-T)-w- (59) 

For the corrections due to real emissions, one has to consider quark-antiquark and two- 
gluon intermediate states, 

jR _ jR _|_ jR 



919 1 9i9 



5.2.1 Quark-antiquark intermediate state 

The starting point here is the quark-antiquark intermediate state contribution to the inclusive 
gluon impact factor (T R = 1/2), 



$ - *° 9q j2^ d ^ T « y 1 ~ —) c~ A m + 



l + ej \C A fc2jfe2 %*k$$ 2 h - h(5 2 f j ' 

(60) 

where /?i and (3 2 are the relative longitudinal momenta (j3i + (3 2 — 1) and fci and fc 2 are the 

— # — * 

transverse momenta {k\ + k 2 — q ) of the produced quark and antiquark, respectively, 
a) quark "inclusive" jet generation 

Taking into account the factor rif arising from the summation over all active quark flavors, 
we have the following contribution: 



i 

t r <*s [ d 2+2e q TV f d ( f ( ( 

a 

V l + eJ\C A ($- kf C (q- k)\q- f ) 2 / ' 

We can split this integral into the sum of singular and non-singular parts. For the singular 
contribution we have 



^ r[i-6] i r 2 (i + 6) (py^-Ur r \\ [d( 



a\ Ca 



x7/,. i I "' s " 



1 + e 



c F r(i + 2 e )r(f- 7 - e )r(n + i + 7 + e ) £ _ 27 



C A T(l + e)r(l - e)r(f - 7 )r(f + 1 + 7 + 2e) 



+ r C 
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Expanding it in e we obtain 
(I R ) 



2ix e(47r) e 



f/.(f) |-/^.(0 



+ c 



-2" 



— + c 27 

C A ^ 



+ n 9 (0 



^x(7,n) + 2C 2 ^h4 



For the regular part of (16"TT) we have 



27r(47r) 



r ( 0i\ C 



xT R [ 1 



2CO C /" ^ 2+2e a 



1 + e/ C J 7T 1+e 
Expanding it in e we get 



2\7~ 



(a 2 ) 



a • I 
n ■ I 



-2-, 



a • n — 1 



fa-n) 2 a - 5 



dC f M CaC 



where we define the function 

h = 7 3 (n,7,C) 



<i 2 a a ■ n — 1 



7T 



a-nr a — 



C 



(62) 



The case of antiquark inclusive generation of the jet is identical to the case of the quark, 
b) both quark and antiquark generate the jet 

The jet momentum is k = k\ + £2 and the jet fraction is 1 = C + C Introducing A as 

h = (k + A , 

the contribution reads 



R a s /pN7- 

m+q 27r(47r) 



k-l) n f f g (a) — 



C A fd 2+2 'A 



d(T R 1- 



2CC 



x 



C F fc 2 







k 2 (Ck + A) ■ [Ck- A 



CU(Cfc + A) 2 (Cfc-A) 2 (£fc + A) 2 (C& - A) 2 A 



(63) 
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In the SCA we need to consider only 

1 A 2 

— m ax — * 



where |A max | = |fc|i2min(£, (). Using again that 

A 2 



f d^A 1 1 (A2 , „ r(l- £ r(l + e) . 2 

7 tt^ A 2 er(l + e) 1 maxJ ~ er(l + 2e) 1 nmx) ' 

9,9+9 2tt e(47r) e T(l + 2e) V / V J NgK ' C F 

l 

J d(T R (mm((X)r (l~Y 



x 





+ e 



a,r(l- e)r (1 + e) - 2 y + -t -y £ 
2tt e 4tt e T l + 2e V / V J NgK J C F 



,23 2, ft 
- - e — + -ln2 
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c) quark "inclusive" jet generation with the antiquark in the jet cone 

— * 

We introduce the vector A such that 
— * — * 

where k coincides with k±, the transverse momentum of the quark generating the jet. 
contribution reads 

i 

jr a * fu2V~^ ft r\ n f rf C-2 7 ^ r f a \ Ca 



(p)^.f)7fc-*»,/. 



m '~ q 2tt(47tY V ) V ~ ) J C ~ JJy \C)C F 

a 

A 2 

_ "max —> 

2CC \ f d 2+2e A 1 



x T R 1 

1 1 + eJ J tt 1 +' A 2 



where now |A max | = \k\R$. We get 

* (py +t -"(k.r) n R 



R _ <*. T(l-e) r 2 (l + e) /r 2 y +e - 2 ir TV r ,2, 
™>-*~~ 2tt e{^Y r(l + 2e) v / v 7 
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2tt e(4 



'MY T(l + 2e){ k ) { ) J C f9 \CjC F 



X( 



C 2l nAP qg {Q 



1 + 2e In i 



Note the overall minus sign, which means that this contribution is a subtractive term to the 
quark "inclusive" jet generation. 



The case of antiquark "inclusive" jet generation with the quark in the jet cone gives the 
same contribution. 



5.2.2 Two-gluon intermediate state 

The starting point here is the gluon-gluon intermediate state contribution to the inclusive 
gluon impact factor, 



+ 



+ 



where (3\ and /3 2 are the relative longitudinal momenta + /3 2 — 1) and fci and k 2 are the 
transverse momenta {ki + k 2 = q) of the two produced gluons. 

a) gluon "inclusive" jet generation 

We need to consider the case of a gluon which generates the jet, while the other is a 
spectator, the case when the other gluon generates the jet being taken into account by a 
factor 2. Thus, we obtain the following integral: 



g;g 



a. 



27r(47r) 
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cHJ c F 
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q- ? 



(q-k) 2 (q~l 



(67) 



The calculation goes along the same lines as in the Section 15.1.11 (case b 2 ). First, we 
separate the C — > 1 singularity, then we add the BFKL subtraction term. Using ()47l) one 
obtains 
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-2 + CC 
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1 1 _L C 2 A: 2 1 2 

(1 - I |) 2 (g - £) 2 + C 2 ((f _ £ )2 ^_ |) 2 j + (1 - £) 2 



r 2+cc 



+ 



+ 



c 2 



<f-f) C (?*-*) 2 (*-§)' 



i C 2 
+ 
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We can split the result into the sum of singular and non-singular parts. For the singular 
contribution we obtain 
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Expanding this result in e we obtain 
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Finally, the e expansion of the divergent part has the form 
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For the regular part, it differs from (1421) only by a factor and reads 
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Expanding it in e one obtains 
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b) both gluons generate the jet 

The jet momentum is k = k\ + &2 and the jet fraction is 1 = C + C- Introducing A a; 

fci = Cfc + A , 

the contribution reads 
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In the SCA we need to consider only 



C (l-C) 



^» d 2 + 2£ A 1 

7T 1+e A 2 



where |A max | = \k\Rmm((X)- Using 
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(71) 



c) gluon "inclusive" jet generation with the other gluon in the jet cone 

— * 

We introduce the vector A such that 

— * — * 

where k coincides with ki, the transverse momentum of the gluon generating the jet. The 
contribution reads 
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where now |A max | = \h\R$. We get 
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Introducing the splitting function P gg , we get 
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5.2.3 Final result for the case of incoming gluon 



Collecting all the contributions calculated in this Section and taking into account the PDFs' 
renormalization counterterm fl28|) and the charge counterterm fl30|) . we find that all singular 
contributions cancel and the result is 



— k 

2vr 



(74) 



y ln 4^ 5(1 - c) 



|p 99 (0 + 2 n/^P^CC) } In ^ - 2C- 27 In P {P gg (C) + 2 n,P w (C)} 

+ CV(1 - |x(n, 7) ^ g + ^ + y + ^ (l + 7 + |) - ^ Q - T) - X 2 (n, 7)) } 

ln(l - 



+ 2C A (1-C 



- ( -2 + K pn C + 



+ C A 

+ 71/ 



X + ^-- 2 + CC 



C (i-0- 



(l + r 27 )x(^,7)-21nC+^/ 2 




C 2 " 



2Cc| + (C 2 + C 2 )(|x(^7) + ^)-^(i-0 



For the Ii,2,3 functions, which enter our final expressions for the quark and gluon contri- 
butions, we obtain the following results: 
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Using the following property of the hypergeometric function, 



2^(1, a, a + 1, C) = a ^ U± [^{n + 1) - ^(a + n) - lnC] C" , 



n=0 



one can easily see that for £ — > 1 



J 2 = 0(lnC), /i = 0(lnC), I 3 = 0(lnC), 



which implies that the integral over ( in ( )56l) and in ( 17^1) is convergent on the upper limit. 



6 Summary 

In this paper we have calculated the NLO vertex (impact factor) for the forward production 
of high-pr jet from an incoming quark or gluon, emitted by a proton, in the "small-cone" 
approximation. This vertex is an ingredient for the calculation of the hard inclusive production 
of a pair of forward high-pr (or Mueller-Navelet) jets in proton collisions. 

At the basis of the calculation of the hard part of the vertex was the definition of NLO 
BFKL parton impact factors; then the collinear factorization (in the MS scheme) with the 
PDFs of the incoming partons was suitably considered. 

We have presented our result for the vertex in the so called {y, n)-representation, which is 
the most convenient one in view of the numerical determination of the cross section for the 
production of a pair of rapidity-separated jets, along the same lines as in Ref. [19J. 

We have explicitly verified that soft and virtual infrared divergences cancel each other, 
whereas the infrared collinear ones are compensated by the PDFs' renormalization countert- 
erms, the remaining ultraviolet divergences being taken care of by the renormalization of the 
QCD coupling. 

In our approach the energy scale sq is an arbitrary parameter, that need not be fixed at 
any definite scale. The dependence on s will disappear in the next-to-leading logarithmic 
approximation in any physical cross section in which jet vertices are used. Indeed, our result 
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for the NLO jet vertex, given by Eqs. (13TT) . ( 132|) . ( 156]) and (!74l) contains contributions ~ ln(s ) 
and these terms are proportional to the LO quark and gluon jet vertices multiplied by the 
BFKL kernel eigenvalue x{ n i v \ This fact guarantees the independence of the jet cross section 
on So within the next-to-leading logarithmic approximation. However, the dependence on this 
energy scale will survive in terms beyond this approximation and will provide a parameter to 
be optimized with the method adopted in Refs. [Tj9]. 

The small-cone approximation, which we adopted here, allows us to obtain explicit analyt- 
ical result for the jet impact factor. In the general case the dependence of the partonic cross 
section on the jet cone parameter has, in the limit R — > 0, the form da ~ A In R + B + 0(R 2 ) 
(see, for instance, [T2] and Appendix C there). In fact, in our work we calculated the coef- 
ficients A and B, neglecting all pieces 0(R 2 ). This can be seen directly from our formulas; 
for example in proceeding from Eq. ( |4"9]) to Eq. f )50|) the contributions O (A 2 ) ~ 0(R 2 ) were 
neglected. The quality of the small-cone approximation has been checked by comparison with 
the results of Monte Carlo calculations which treat the cone size exactly, both for the cases 
of unpolarized and polarized jet cross sections. Very good agreement between the results of 
the small-cone approximation and the Monte Carlo calculations was found even for cone sizes 
of up to R = 0.7, see [TB] for more details and references. Therefore having this experience 
with the jet production in Bjorken kinematics s ~ Q 2 , there is the hope that the small-cone 
approximation could also be an adequate tool for describing Mueller-Navelet jets for an exper- 
imentally relevant choice of the jet cone parameter, R ~ 0.5. Another important application 
of the small-cone approximation method could be the possibility to perform a semi-analytical 
check of the complicated numerical approaches to Mueller-Navelet jet production which treat 
the cone size exactly, in the limit of small values of R. 
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A Appendix 



We list here some useful integrals 
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In the integrals below, / 2 = is assumed 
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